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ABSTRACT 

Context. We have analysed far-infrared spectra of 32 circumstellar disks around Herbig Ae/Be and T Tauri stars obtained within 
the Herschel key programme Dust, Ice and Gas in Time (DIGIT). The spectra were taken with the Photodetector Array Camera and 
Spectrometer (PACS) on board the Herschel Space Observatory. In this paper we focus on the detection and analysis of the 69 pm 
emission band of the crystalline silicate forsterite. 

Aims. This work aims at providing an overview of the 69 pm forsterite bands present in the DIGIT sample. We use characteristics 
of the emission band (peak position and FWHM) to derive the dust temperature and to constrain the iron content of the crystalline 
silicates. With this information, constraints can be placed on the spatial distribution of the forsterite in the disk and the formation 
history of the crystalline grains. 

Methods. The 69 pm forsterite emission feature is analysed in terms of position and shape to derive the temperature and composition 
of the dust by comparison to laboratory spectra of that band. The PACS spectra are combined with existing Spitzer IRS spectra and 
we compare the presence and strength of the 69 pm band to the forsterite bands at shorter wavelengths. 

Results. A total of 32 disk sources have been observed. Out of these 32, 8 sources show a 69 pm emission feature that can be attributed 
to forsterite. With the exception of the T Tauri star AS 205, all of the detections are for disks associated with Herbig Ae/Be stars. 
Most of the forsterite grains that give rise to the 69 pm bands are found to be warm (~ 100-200 K) and iron-poor (less than ~2% iron). 
AB Aur is the only source where the emission cannot be fitted with iron-free forsterite requiring approximately 3-4% of iron. 
Conclusions. Our findings support the hypothesis that the forsterite grains form through an equilibrium condensation process at high 
temperatures. The large width of the emission band in some sources may indicate the presence of forsterite reservoirs at different 
temperatures. The connection between the strength of the 69 and 33 pm bands shows that at least part of the emission in these two 
bands originates fom the same dust grains. We further find that any model that can explain the PACS and the Spitzer IRS observations 
must take the effects of a wavelength dependent optical depth into account. We find weak indications of a correlation of the detection 
rate of the 69 pm band with the spectral type of the host stars in our sample. However, the sample size is too small to obtain a definitive 
result. 

Key words, star: Herbig Ae/Be stars - star: T Tauri stars - forsterite - infrared - young stellar objects - spectroscopy 



1. Introduction 

The diversity of planetary system architectures and planet prop- 
erties is certainly related to a range of physical and chemica l 
conditions in protoplanetary disks (e.g Morda sini et al.ll2012h . 
Dust particles in disks undergo a wide range of physical and 
chemical processes, ranging from dust growth through cohesive 
coagulation to annealing, sublimation and re-condensation (e.g . 
Natta "etat]|2007t iHenning & Meeu3l201ll: iGail & Hoppi201(| 
Gail 2010). Infrared (IR) spectroscopy turns out to be a versa- 
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tile tool to characterise the physical and chemical structure of 
the dust particles through the study of emission features arising 
from small dust grains. 

Near- and mid-IR radiation comes from the upper optically 
thin layer of the disk atmosphere heated by stellar radiatio n 
dMen'shchikov & Henningl 119971: IChiang & Goldreichl 1 19971) . 
The infrared spectra of these regions i n the disks have been 
intensively studied , first by ISO (e.g. iBouwman et al.l [2003: 
Meeus et al.l 120011) for intermediate-mass Herbig Ae/Be stars 
and later with much improved sensitivity provide d by Spitzer for 
disks around brown dwarfs/low-mass stars (e.g. iPascucci et al.l 
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20091) T Tauri stars (IBouwman et alJl2008[lKessler-Silacci et al l 
2006; lOlofsson et al.l 120091: iFurlan et al.l l2009t IWatson etafl 
20091) and Herbig Ae/Be stars (lJuhasz et al.l l2010h . Regions 
closer to the mid-plane can only be probed by observations at 
(sub)-millimetre wavelengths where the disks become optically 
thin. 

These deeper layers of protoplanetary disks are optically 
thick to instruments in the near and mid-IR range. At around 
10 pm, the emission of regions closer to the m idplane is de- 
tectab le even in relatively massive disks (see, e.g. lMulders et alj 
An overview of the spectral feature s at these wavelength s 
and the information they carry is given bv Ivan Dishoeckl d2004). 
Studie s in this regime have been conducted wit h ISO, for exam- 
plebv lMalfait et all dl998l) . lMalfaitl d 1999b and lLorenzetti et all 
(120021) . 

The emission band at -69 pm, associated with crystalline 
olivine grains, is much more sensitive to the temperature and 
the iron content of the crystalline dust parti cles than any of the 
olivine bands at shorter waveleng ths (e.g. iBowev et~al1 f2002: 
iKoike et alJl2003HSuto et alJ 2006). Prior to the Herschel obser- 
vations, the 69 pm band has only been detected by ISO in one 
protoplanetary disk (HD 100546, Mal fait et al.ll9 98). However, 
that detection had a low signal-to-noise ratio (S/N) and the 
shape of the emission band was never analysed in terms of tem- 
perature or iron content. Herschel PACS spec tra of that particular 
source were discussed bv lSturm et al.l (12010). who presented the 
69 pm band at unprecedented S/N and resolution and who found 
the emission to emerge from very iron-poor dust at a temperature 
of about -150-200 K. 

In this paper we present an overview and analysis on the 
69/im forsterite emission band in the disk sample observed 
within the Herschel key programme DIGIT. The sources in the 
DIGIT sample cover a wide range of disk properties including 
age, luminosity and effective temperature. Through the chosen 
set of disk properties the sample includes examples of different 
evolutionary phases. All sources in the brightness-limited sam- 
ple have been drawn from previous studies with a focus on those 
objects for which high-quality Spitzer IRS 5-35 pm spectra ex- 
ist. This allows for a comparison of the hot, geometrically thin 
surface layer to the interior of the disks. All Herschel observa- 
tions in this work were taken with the PACS instrument. 

In this paper we search for the 69 /mi emission band of 
forsterite (see section[3). In section|4]we analyse the position and 
shape of the detected emission bands in terms of temperature and 
iron content. A brief discussion on the possible influence of grain 
size is included. Furthermore, we compare the 69 pm band to 
the forsterite emission bands at mid-IR wavelengths (e.g. 16 and 
33 pm) from Spitzer observations. We search for a relation in the 
peak over continuum ratio in the different wavelength regimes, 
which would give further information about the spatial distribu- 
tion and the formation history of the forsterite in the disk. 



2. Observation and data reduction 

In this work we p resent Photodetector Array Camera and 
Spectrometer ( PAC S : IPo glitsch et al .120 1 01) spectra of protoplan- 
etary disks in a sample of Herbig Ae/Be and T Tauri systems. 
These observations were taken as part of the Dust, Ice and Gas 
in Time (DIGIT) Herschel key programme. The list of sources 
we observed is provided in Table [1] 

The PACS instrument consists of a 5x5 array of 9.4"x9.4" 
spatial pixels (here after referred to as spaxels) covering the 
spectral range from 51-210 pm with A/AA ~ 1000-3000. The 



spectrum is divided into four segments, covering A ~ 50-75, 70- 
105, 100-145, and 140-210 pm. The spatial resolution of PACS 
ranges from ~9" at 50 pm to ~ 18" at 210 pm. All of our targets 
were observed in the standard "rangescan" spectroscopy mode 
with a grating stepsize corre sponding to Nyquist sampling (see 
further iPoglitsch et al.l2010l) . In this paper we focus on the spec- 
tral range of 67-72 pm observed in the second spectral order 
with the "blue" detector. 

We processed our data u sing the Herschel Interactive 
Processing Environment (HIPE, Ott 2010) using calibration ver- 
sion 42 and standard pipeline scripts. The infrared background 
emission was removed using two chop-nod positions 6' from 
the source in opposite directions. Absolute flux calibration was 
achieved by normalising our spectra to the emission from the 
telescope mirror itself as measured by the off-source positions, 
and a detailed model of the telescope emission available in 
HIPE. The sources were usually well-centred on the central 
spaxel which we used to extract the spectra as this provided the 
highest signal to noise ratio (S/N) in the spectra. However, small 
pointing errors and drifts of the telescope can lead to flux losses 
and spectral artifacts. To mitigate this we scaled the spectra de- 
rived from the central spaxel to the integrated flux over the entire 
array. This approach guarantees the best absolute flux calibration 
with the highest S/N spectra. Only for one case, HD 142666 , the 
telescope had such a large mispointing that the central spaxel did 
not contain most of the source flux. Here we opted to use the 
spectra integrated over the central 3x3 spaxels to recover the to- 
tal flux. Spectral rebinning was done with an oversampling of a 
factor of two and an upscaling of a factor of one corresponding 
to Nyquist sampling. For further details on the data reduction 
procedure for the entire DIGIT dataset we refer to Green et al . 
(subm.). 



3. Detection of forsterite 

3.1. Characteristics of the 69 pm emission band 

Olivine, a nesosilicate with an orthorombic crystal structure, 
forms a complete solid solution series from forsterite (Mg 2 Si04) 
to fayalite (Fe2Si04). The general chemical composition is given 
by Mg2(i- A )Fe2 A Si04, where x is the fraction of iron cations 
relative to magnesium. Magnesium-rich olivines give rise to a 
characteristic emission band at ~69 pm. The most pronounced 
feature is found in forste rite -the iron-free end member of the 
olivines dHenningl uOlO). The position and shape of the band 
depends on, in order of importance, the iron content, the temper- 
ature and the shape and size of the dust grains. 

In this paper we will use the term "forsterite" or 
"pure forsterite" for completely iron-free olivines. "Iron-poor 
forsterite" is used for olivines with less than 5% of iron cations 
relative to magnesium. When using the term "iron-poor olivines" 
we refer to olivines with no more than 10% of iron. These terms 
are used to refer to the dust particles; the emission band is al- 
ways labelled "forsterite emission band" or "69 pm band" as it 
is strongest in pure forsterite and vanishes completely in iron- 
rich olivines. 

3.1 .1 . Effects of iron content 

For olivines in general, the position and shape of the emission 
peak strongly depends on the iron content. The 69 pm band is 
shifted linearly toward longer wavelen gths with the amount of 
iron in the crystals (e.g. Koike et al. 2003). Even a small fraction 
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Coordinates (J 2000) 



Star 


RA [h m s] 


DEC [d m s] 


spectral type 


ref. 


log(LJL ) 


ref 


age (Myr) 


ref 


HD 150193 


16 40 17.92 


-23 53 45.2 


Al-3Ve 


(1) 


1.47 


(11) 


>2.0 


(11) 


HD 97048 


11 08 03.32 


-77 39 17.5 


B9.5Ve+sh 


(1) 


1.61 


(11) 


>2.0 


(11) 


HD 169142 


18 24 29.78 


-29 46 49.4 


A5V 


(5) 


1.58 


(5) 


7.7+2.0 


(10) 


HD 98922 


11 22 31.67 


-53 22 11.5 


B9 Ve 


(11) 


>2.96 


(11) 


<0.01 


(15) 


HD 100453 


11 33 05.58 


-54 19 28.5 


A9 Ve 


(1) 


0.95 


(12) 


>10 


(10) 


HD 135344 


15 15 48.95 


-37 08 56.1 


F4Ve 


(1) 


1.02 


(7) 


9+2 


(7) 


HD 179218 


19 11 11.25 


+ 15 47 15.6 


B9e 


(1) 


2.50 


(11) 


1.1+0.7 


(14) 


HD 203024 


21 16 03.02 


+68 54 52.1 


A/B8.5 V 


(19) 


-2.0 


(19) 






IRS 48 


16 27 37.19 


-24 30 35.0 


A? 


(2) 


1.16 


(8) 


-15 


(8) 


SR21 


16 27 10.28 


-24 19 12.5 


G2.5 


(3) 


1.45 


(3) 


-1.0 


(3) 


HD 38120 


5 43 11.89 


.4 59 49.9 


Be* 


(2) 


1.74 


(9) 


<1.0 


(13) 


HT Lup 


15 45 12.87 


-34 17 30.6 


K3 Ve 


(2) 






0.5+0.1 


(16) 


HD 35187 


5 24 01.17 


+24 57 37.6 


A2Ve+A7Ve 


(1) 


1.24 


(10) 


9.0±2 


(10) 


S CrA 


19 01 08.60 


-36 57 20.0 


K6e 


(3) 


0.36 


(3) 


-3.0 


(3) 


HD 104237 


12 00 05.08 


-78 11 34.6 


A7.5Ve + K3 


(4) 


1.55 


(4) 


5.5+0.5 


(10) 


RYLup 


15 59 28.39 


-40 21 51.2 


GOV 


(2) 






29.7+7.5 


(16) 


HD 144432 


16 06 57.96 


-27 43 09.8 


A9IVev 


(1) 


>1.48 


(11) 


5.5+2.0 


(14) 


AS 205 


16 11 31.35 


-18 38 26.1 


K5 


(3) 


0.85 


(3) 


-0.1 


(3) 


HD 144668 


16 08 34.29 


-39 06 18.3 


A5-7III/IVe+sh 


(1) 


1.71 


(10) 


2.8+1.0 


(10) 


RULup 


15 56 42.31 


-37 49 15.5 


G5 Ve 


(2) 






55.3+10.4 


(16) 


HD 141569 


15 49 57.75 


-3 55 16.4 


AO Ve 


(1) 


1.35 


(11) 


4.7+0.3 


(10) 


EC 82 


18 29 56.89 


+ 1 14 46.5 


K8D 


(18) 


-0.179 


(18) 


-1.5 


(18) 


HD 50138 


6 51 33.40 


-6 57 59.4 


B9 


(2) 


2.85 


(11) 


0.5+0.2 


(16) 


HD 142666 


15 56 40.02 


-22 01 40.0 


A8Ve 


(1) 


1.13 


(10) 


9.0+2.0 


(10) 


DGTau 


4 27 04.70 


+26 06 16.3 


G Ve 


(2) 


-0.05 


(17) 


-0.6 


(17) 


RNO 90 


16 34 09.17 


-15 48 16.8 


G5D 


(2) 










HD 100546 


11 33 25.44 


-70 11 41.2 


B9Vne 


(1) 


1.51 


(11) 


>10 


(11) 


AB Aur 


4 55 45.84 


+30 33 04.3 


AOVe+sh 


(1) 


1.68 


(11) 


5.0+1.0 


(10) 


HD 163296 


17 56 21.29 


-21 57 21.9 


A3Ve 


(1) 


1.48 


(11) 


5.5+0.5 


(10) 


HD 142527 


15 56 41.89 


-42 19 23.3 


F7IIIe 


(1) 


1.84 


(11) 


2.0+0.5 


(10) 


HD 36112 


5 30 27.53 


+25 19 57.1 


A5 Ve 


(2) 


1.35 


(11) 


3.7+2.0 


(10) 


HD 139614 


15 40 46.38 


-42 29 53.5 


A7Ve 


(1) 


0.88 


(10) 


10+2.0 


(10) 



The references fo r spe ctral class ifications and a ges are: (l )lAck e & van den AnckeH d2004h . (2 ) The Simbad datab ase, ( 3)lPrato et al.ld2003h. (4 ) 
Boh m et al.l d2004h. (5)Pfilling et all (120121). (6)lMeeus et all j2002h. (7)lMuller et al.l (hoi lh . (8)lBrown et al.l d2012l) . ( 9)lHernandez et alj J2 005). 
(TO)lMeeus et alj d2012n. (Tl)lvan den Ancker et all d!998h. Vl2)lDominik et alJ d2003h. ("l3)ISartori et all <|2010h. a4)lFolsom et al] d2012h . (15) 
Man oi et al.1 d2006h . (16) iTetzlaff et alj d201 ll) . (17) lPalla & StahleJ ( 120021) . (18) lwinston et alj d2010l) . q9) lMiroshnichenko et alj d!999l) . 
Note: Luminosities and ages are derived by different methods. 



of iron (x - 0.05) leads to a band position at A > 70 fim (see Fig. 
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Fig. 1. The influence of the iron cation fraction relative to mag- 
nesium in the olivine material on the peak wavel ength of the 
69 um band. The data (black crosses) are from iKoike et all 
(2003), taken at 295 K and the blue line is a linear fit. 



Laboratory measurements are only available for pure 
forsterite and for olivines with more than ~10 % of iron con- 
tent. From the available data (presented in Fig. [TJ a linear con- 
nection between iron content in the olivine dust grains and the 
peak position of the 69 /urn band can be derived. Most of the 
available laboratory measurements were taken at room tempera- 
ture. However, corrections for temperature effects can be applied 
if necessary, based on the temperature dependence of the pure 
forsterite and a sample with 9.3 % iron (Koike et al. 2006). 

The FWHM of the emission band is also increasing with the 
iron content in the grains. However, the available data is sparse, 
introducing additional uncertainties in the quantitative analysis 
of our detections. Therefore our estimated upper limits on the 
iron content in the crystalline silicates from which the 69 /im 
band emerges are mostly based on the position of the peak. 

3.1.2. Effects of temperature 

Next to the strong dependency on the iron content, the posi- 
tion and width of the 69 fim forsterite band also depends on the 
grain temperature. This temperature dependence h as been char- 
acterised in several laboratory experim ents (e.g. iBowev et al.l 
l2002HKoike et alJ2006r.lSuto et alJl2006h . covering a range from 
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8 to 295 K. In Fig. [2] we show the FWHM and peak position of 
the 69 //m band at different temperature s. The mass abso rption 
coefficients from the optical constants bv lSuto et al.l (12006) were 
computed u sing "distributio n of hollow spheres" scattering the- 
ory (DHS, IMin et al.ll2003h . We used a filling factor / max of 1 
and a grain size of 0. 1 /mi. From 8 to 295 K the peak position 
moved from 68.8 /mi to 69.6 //m and the FWHM of the emis- 
sion band increases from ~0.2 to 1.1 //m as can be seen from 
Fig. [2] Note that all of these measurements were carried out on 
pure forsterite. The temperature dependence of the 6 9 /mi band 
at 9.3% iron content is shown in Koi ke et a n j2"006h . The effect 
is qualitatively similar to that in pure forsterite. 



(lJuhasz et al.l20 10). The temperature range covered is 5 0-295 K 
and the data is based on the results bv lSuto et al.l (120061) . 
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Fig. 2. The effect of temperature on the peak position and 
FWHM of the 69 ji va band. The blue crosses are based on op- 
tical constants from Sut o et al. I (I2006I) using the distribution of 
hollow spheres bv lMin et al.l d2003h with a grain size of 0.1 /mi. 
The green asterisk s represent mass absorption coefficients from 
Koik e et alJ ((2006). Both datasets were taken on pure forsterite. 



The offset between the different datasets shown in Fig. [2] 
can be explained by the effects of grain size, grain shape, 
amount of lattice distor tions, possible effec ts from the scatter- 
ing theory (see e.g., iMutschke et al.l 120091) as well as even- 
tual environmental contamin ation during laboratory experiments 
dHenn ing & M utschkdEo iO). Such differences can be the effect 
of different methods used in sample preparation. In the follow- 
ing section we will discuss, and take into account, the effects of 
gr ain size and shape b ut not that of lattice distortions. As noted 
by iKoike et aD d2010h . the number of lattice distortions affects 
the shape (and to a minor extent the position) of the emission 
band. The broade ning of the IR emission bands observed by 
Koike et al. ( 2010) and attributed to lattice distortions, however, 
is too small to be detected in the presence of the other effects 
we discuss in this paper. Therefore, we will will not take it into 
account in our further analysis. 

3.1 .3. Effects of grain shape 

Another factor that has to be taken into account is the shape 
model which is needed to compute the mass absorption coef- 
ficients of dust grains from t he refractive index of the material 
dHenning & Mutschke1 l20Tol) . In Fig. [3] we compare the calcu- 
lated peak position and FWHM for compact spherical grains 
(Mie scattering ) and the distribution of hollow spheres (DHS; 
IMin et all 2003) model with a filling factor / max = 1, which was 
found to be a good model for the crystalline silicate grains pro- 
ducing features in the Spitzer IRS wavelength regime (5-35 /mi) 
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Fig. 3. The effect of the shape model used for the forsterite grains 
on the peak position and FWH M. Shown as bl ue crosses is the 
distribution of hollow spheres (IMin et al. 2003) while green as- 
terisks stand for Mie scattering on spherical grain s. Both exam- 
ples are based on the data from lSuto et al.l (I2006I) and cover the 
temperature range from 50 to 295 K. The grain size is 0. 1 /mi 
and the black arrow indicates the temperature dependency. 

On average, both the peak position and the FWHM com- 
puted with the DHS model are shifted to larger values by about 
0.1/mi compared to Mie scattering. However, as can be seen 
from Fig. [3] the grain shape has little influence on the relation 
between the grain temperature and the peak position and FWHM 
of the forsterite band, which shows the far stronger effect. The 
same holds for the dependency on the iron content. From mid-IR 
observations we know that purely spherical grains do not repro- 
duce the observed spectral features well. In this work, therefore, 
we adopt the DHS model with a maximum filling factor / max = 1 
as shown in Fig. [3] Results from a third wi dely used theor y, the 
continuous distribution of ellipsoids (CDE; IMin et aLll2003l) are 
almost identical to those obtained from the DHS model. The dif- 
ference is too small to lead to significant changes in the results 
for temperature or iron content. Also, as no grain size effects can 
be taken into account in this model, we do not consider CDE in 
the remainder of this paper. 

3.1 .4. Effects of grain size 

A final parameter we discuss here is the size of the dust particles 
which also has an effect on the shape and peak position of the 
69 /mi band. Grain size effects could be important as we can ex- 
pect to observe dust that is deeply embedded in the disk where 
grain growth and settling will likely play a role. We computed 
mass absorption coefficients for 7 different grain sizes, between 
1 and 50/m i, using the DHS scattering theory and optical con- 
stants from ISuto et ail (|2()()6). and the results are displayed in 
Fig. 4. The changes in the band can be clearly seen for grains 
larger than 1 /mi, while the band profile is almost indistinguish- 
able for grains with radii between 0. 1 and 1 .0 //m. The maximum 
shift in the peak position we expect due to changes in the grain 
size is in the order of 0.2 /mi. The effects seen in the FWHM 
of the forsterite band is larger, however, the peak over contin- 
uum ratio decreases significantly for grains larger than 10 /mi. 
With the S/N we achieve in our sample it would be impossible 
to detect the emission of grains larger than about 15-20 //m in ra- 
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Fig. 4. The effect of grain size on band shape. The top panel 
shows the profile of the 69 /mi forsterite band at 150 K for dif- 
ferent grain siz es. The profiles w ere computed using the optical 
constants from Suto et al. (2006) and the distribution of hollow 
spheres ( Min et al. 2003). The lower panel shows the relation be- 
tween peak position and FWHM of the profiles shown in the top 
panel. The asterisks in the lower panel correspond to the grain 
sizes in the upper panel, starting with 1 /mi at the bottom. 

dius, so we do not expect to detect the very broad emission bands 
which can be associated with large grains. The typical grain sizes 
we consider in this study are 0. 1-1 .0 /mi. However, in our analy- 
sis (see section l4~4l we test the conformity of the detected bands 
with larger grains. 

3.1 .5. Summary of effects 

In summary, the iron content and the grain temperature have the 
largest effect on the 69 //m feature. However, for small measured 
shifts of the order of 0.2 /zm in the peak position and FWHM, a 
degeneracy between the parameters exists and no clear distinc- 
tion can be made between small changes in iron content, temper- 
ature, or grain size and shape. In our analysis we will discuss the 
constraints on iron content and grain temperature, and discuss if 
we see any indications for substantial grain growth beyond the 
typical grain size of 0. 1-1 /mi, the grain sizes typically observed 
in the mid-IR, taking into account these possible degeneracies. 

3.2. Searching for the 69 //m band in the PACS spectra 

To search for and analyse the 69 /mi forsterite band in our sam- 
ple, we examine the 67-71. 5 pm region in the spectra, shown 
in the left column of Figures IA.1I to IA.8I in the Appendix. 
Aside from a few sources with very strong emission bands (e.g. 



HD 100546) noise reduction is important to separate forsterite 
bands from narrow features such as gas lines or noise peaks. 

In order to reduce the noise and to concentrate on wide fea- 
tures, all spectra underwent noise-filtering in their 67-7 1.5 /mi 
range. We compute the Fourier-transform (FT) of each spectrum 
and remove the contribution of all frequencies outside of the low 
frequency peak. After inverse FT the noise-filtered spectrum is 
overplotted on the unmodified data (see left column of Figures 
IA.H4A.8l l. In general the filtering greatly improved the visibility 
of wider features such as the forsterite bands while gas lines, 
noise and other narrow spikes in the spectra are suppressed. 
However, some very strong gas lines (e.g. the 71 /mi OH fea- 
ture in DG Tau) do not completely vanish through the filtering . 
Therefore, we carefully compared all forsterite detections with 
the unmodified spectra to avoid confusion with emission bands 
other than the forsterite 69 /mi feature. 

We use the noise-filtered spectra to search for possible 
69 //m forsterite emission bands. In some cases the noise-filtered 
spectra suffer from the Gibbs-phenomenon at the outer edges, 
introducing artificial features that could be confused with emis- 
sion bands. To avoid complications we restrict the position for 
the peak of possible emission bands to the 68-7 1/mi range. 
Visual inspection did not reveal any forsterite peaks outside this 
range. Some of the spectral resolution is lost due to the noise- 
filtering process. To make sure that no false positives were cre- 
ated through the filtering we compare the filtered to the unmod- 
ified spectrum. The comparison is done by overplotting the un- 
modified with the filtered version of the spectra (see leftmost 
column in Figures lA~TT4A.8t . 

The 67-7 1 .5 //m region of each noise-filtered spectrum is fit- 
ted with a 2 nd or 3 rd order polynomial to account for the contin- 
uum emission, and a superimposed Lorentz profile to describe 
one possible emission band. The Lorentz profile is restricted to 
a peak position in the range of 68-71 /mi, a FWHM of less than 
1 .5 /mi. We considered only positive values for the integrated 
flux to avoid fitting absorption bands. These restrictions prevent 
the fit from converging on gas lines or features that are extended 
outside the fitted spectral range. Any fit solution where one or 
more parameters reached these limits was not used to claim a 
detection. 

As the data reduction does not provide reliable uncertainties, 
the first fit is done using equal weights (1) for all points. This 
leads to improbably low reduced^- 2 (y 2 ) values, so the fit is re- 
peated with assumed error bars of y^ 2 for each datapoint. If 
the parameter values remain the same, this will result in a new 
X 2 r = 1 • If the parameters change, the procedure is iterated until 
0.95 < xl < 1 -05 is reached. A proper estimate of the uncertain- 
ties in the spectrum is important to obtain reliable estimates for 
the uncertainties of the model parameters. 

The best-fit models overplotted on the noise-filtered spectra 
are shown in the middle column of Figures lA~TT4A.8l The param- 
eter values found for the Lorentzian are listed in Table [2] The 
table also includes the formal uncertainties of the peak position 
and the FWHM derived by the IDL implementation of the least 
squares minimisation mpfit ( Markwa rdtl2009l) from the covari- 
ance matrix. The uncertainty of the integrated flux is computed 
from the standard deviation of the residuals. 

Finally we check how likely our best-fit Lorentz curve is to 
be explained by a combination of residual noise. The best-fit 
model is subtracted from the noise-filtered data. We then take 
the standard deviation of the residuals as an estimate for error 
bars (in which case the residuals could be fitted with F(A) = 0) 
and determine the probability that the best-fit Lorentzian is a re- 
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suit of that noise. We compute 4. Analysis of forsterite emission bands 

„ 4.1. The iron content and temperature of the crystalline 

D = )^ (l/icr^m, P )f (1) silicates 



where i is the index of the wavelength points in the spectrum, L 
the Lorentz function and p the parameters from the fit. 

Assuming D follows ax 2 distribution with N (= number of 
wavelength points minus number of parameters) degrees of free- 
dom, we compute the probability P to find a value of D as high 
as measured or higher. The values of P are listed in Table[2] The 
residuals, overplotted with the F(A) = and the error bars as well 
as the best-fit Lorentzian is shown in the right column of Figures 

The criteria for a detection of the forsterite 69 yum band are aj 
integrated flux/error ratio > 3 in the best-fit model, a low prob- 
ability for the fitted band to be a result of residual noise (we 
do not have to specify a precise number as the probability is in 
all cases either almost unity or almost zero), and that none of 
the fit parameters have reached their boundaries. Furthermore, 
all band models that meet these criteria are checked against the 
unmodified spectrum to make sure that the band was not sig- 
nificantly altered through the noise-filtering. These checks were 
done through over-plotting the best-fit model on the unmodified 
spectrum. A more rigorous solution would be to fit the model 
also to the unmodified spectrum and check if the parameter val- 
ues are the same as with the filtered spectrum. As the emission 
band is very weak in some cases (e.g. AS 205, HD 104237) this 
would require some fine tuning of start values for the parame- 
ters. Therefore, we only visually checked the model overplotted 
on the unmodified spectrum. 

The sources with firm detections following the analysis pre- 
sented in this paper are AB Aur, HD 100546, HD 104237, 
HD 141569, HD 179218, HD 144668, AS 205 and IRS 48. 
Continuum-subtracted 67-72yum spectra of these objects are 
shown in Figure|8] The detections in AS 205 and HD 144668 are 
weak and could be labelled as marginal. However, all the formal 
criteria described above are fulfilled and the visual inspection of 
the spectrum in the 67-71 pm range (see Figures lA~TT4A.8l in the 
appendix) shows that the fitted peak is stronger than any other 
feature in its vicinity. Several other sources also show signs of 
an emission band in the region of 69 pm but did not allow for a 
definitive identification. We will discuss these cases below. 

The upper limits for the sources without detectable bands are 
computed by inserting a band of the same shape and position as 
the one found in HD 100546 and scale the integrated flux until 
3-3.2 <x over the residuals is reached. 

3.3. False positives 

For six systems, namely DG Tau, HD 100453, HD 203024, 
HD 35187, HT Lup and SR 21, our Lorentz curve fitting for- 
mally fits a peak in the studied wavelength range. However, the 
cr value of the fits is low and the probability to find a value in a 
X 2 distribution greater than or equal to the one found for com- 
paring the Lorentzian to the noise in the residuals is for these 
six systems (near) 100 %. This means that the probability for the 
central peak to be the result of residual noise is very high and we, 
therefore, classify these as false positives. By visually inspecting 
the spectra of these sources, one can clearly see multiple "peaks" 
not associated with forsterite emission and which are most likely 
low level spectral artifacts. In the following analysis we will treat 
the values quoted in Table|2]for these six sources as upper limits. 
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Fig. 5. FWHM of the forsterite '69 pm band' plotted over the 
peak position, fitted with a Lorentz profile. The red squares 
with error bars represent measurements taken from our PACS 
observations. The numbers refer to the stars as follows: (0) 
HD 100546, (1) HD 104237, (2) ABAur, (3) HD 141569, (4) 
HD 179218, (5) IRS 48, (6) HP 144668, (7 ) AS 205. The green 
asterisks are taken fromlKoike et al.l (120031) and the blue crosses 
are based on ISuto et al.l (120061) . computed with the DHS model 
dMin et al.ll2003l) and a grain size of 0.1 fim. The black arrow 
indicates the temperature dependency. 

As discussed in Sections 13.1.11 and 13.1.21 the two parame- 
ters most strongly influencing the forsterite band are the iron- 
to-magnesium ratio and the grain temperature. In Figure |5]we 
plot the measured peak position and FWHM of the detected 
69 yum features. Comparing these observed values to the labo- 
ratory measurements shown in Figures [T]and|2]it is immediately 
clear that the crystalline olivine grains can not contain much 
more than a few percent iron and have a typical temperature of 
about 150 K. Both AS 205 and AB Aur show deviating values 
compared to the other sources. Where the observed 69 pm band 
of AS 205 seems to be consistent only with pure forsterite grains, 
the band position observed for AB Aur can only be explained if 
the olivine grains contain some fraction of iron. 

Ideally, one could have used extensive laboratory studies, ex- 
ploring in detail the combined effects of iron contend and grain 
temperature on the forsterite band. Unfortunately, almost all of 
the laboratory measurements of the position of the forsterite 
band as a function of ir on content were taken at room temper- 
ature dKoike et alj |2003). The only temperature-dependent labo- 
ratory measure ments available are for pure forsterite (M gaSiO/Q 
crystals (e.g. ISuto et al.l l2006t IKoike et alJ 120061) and for 
iron-poor olivine with a relative iron mass fraction of ~10% 
(Koike et al. 2006), for which the peak of the emission band has 
already shifted to ~72.2 pm, far beyond any of the observed peak 
positions. 

To make a quantitative statement about the iron content and 
temperature of the olivine grains, we assume that the peak posi- 
tion and FWHM of the 69 fim feature scale linearly with the two 
main parameters. These assumptions are consistent with the data 
shown in Figures[T]and|2] We fitted Lorenzians to the l aboratory 
measurements of olivine bands by IKoike etaO d2003l) and lin- 
early interpolated between those, creating a fine grid of feature 
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Star 


Peak 


FWHM 


Flux 




probability 




[pm] 


[jim] 


10-' 6 [W/m 2 ] 




for* 2 [%] 


Detections 


AB Aur 


69.957 


±0.018 


0.509 ± 0.063 


10.14 ±0.90 


11.2 





HD 100546 


69.194 


± 0.004 


0.681 ± 0.015 


94.60 ± 1.27 


74.6 





HD 104237 


69.224 


± 0.012 


0.351 ± 0.038 


2.14 ± 0.21 


9.8 





HD 141569 


69.303 


± 0.021 


0.600 ± 0.092 


6.32 ± 0.49 


12.9 





HD 179218 


69.196 


± 0.005 


0.502 ± 0.020 


7.33 ±0.18 


39.7 





HD 144668 


69.088 


± 0.014 


0.599 ± 0.057 


2.60 ±0.13 


19.8 





IRS 48 


69.168 


± 0.007 


0.521 ± 0.025 


6.11 ± 0.21 


28.9 





AS 205 


68.745 


± 0.017 


0.595 ± 0.085 


4.59 ± 0.29 


15.6 





False positives 


DGTau 


69.187 


± 0.038 


0.506 ± 0.153 


2.21 ±0.39 


5.6 


100 


HD 100453 


69.427 


± 0.043 


0.749 ± 0.157 


2.93 ± 0.34 


8.6 


100 


HD 203024 


69.466 


± 0.043 


0.712 ± 0.204 


3.26 ± 0.41 


8.0 


100 


HD 35187 


68.605 


± 0.039 


0.695 ±0.153 


1.98 ± 0.23 


8.5 


100 


HT Lup 


69.723 


± 0.016 


0.285 ± 0.056 


3.33 ± 0.61 


5.5 


100 


SR21 


69.390 


± 0.025 


0.510 ±0.101 


2.60 ±0.31 


8.4 


94 


Non-detections (upper limits) 


HD 144432 








0.84 ± 0.27 


3.1 


100 


HD 97048 








1.96 ± 0.63 


3.1 


100 


HD 36112 








0.72 ± 0.23 


3.1 


100 


HD 142666 








6.53 ±2.11 


3.1 


100 


HD 163296 








0.82 ± 0.26 


3.1 


100 


HD 50138 








2.08 ± 0.67 


3.1 


100 


HD 150193 








0.67 ± 0.22 


3.1 


100 


HD 135344B 








1.12 ± 0.36 


3.1 


100 


HD 169142 








1.73 ± 0.56 


3.1 


100 


SCrA 








1.21 ± 0.39 


3.1 


100 


HD 139614 








0.75 ± 0.24 


3.1 


100 


HD 98922 








1.15 ±0.37 


3.1 


100 


RNO 90 








1.06 ± 0.34 


3.1 


100 


RULup 








0.90 ± 0.29 


3.1 


100 


RYLup 








0.98 ± 0.32 


3.1 


100 


EC 82 








5.09 ± 1.64 


3.1 


100 


HD 38120 








0.60 ±0.19 


3.1 


100 


HD 142527 








2.67 ± 0.86 


3.1 


100 



Description: Fitted-band centre position, FWHM, flux and flux/uncertainty as well as the probability of finding a value in a* 2 distribution greater 
than or equal to the one found for comparing the Lorentzian to the noise in the residuals (see description in the text). The limits for the parameters 
of the Lorentzian were as follows: 68 yum < Peak < 71 yum, FWHM < 1.5 /jm and 10~ 30 W/m 2 < flux. The first part of the table lists the firm 
detections. In the second part all sources with an emission band that has a formal significance of more than three cr but are rejected for various 
reasons (see text for a discussion) are shown. The third part contains all sources where no band could be fitted and we present upper limits instead. 



shapes as a function of iron fraction (between 0-6 %). We also 
fitted the position and FWHM o f the 69 /j.m band as function of 
grain temperature to the data by Koike et al. (2006); Suto et al] 
(2006) in the range of ~ 10-300 K and linearly interpolated be- 
tween them. 

The linear combination of both the iron and the temperature 
dependence provides us with a model grid for the band charac- 
teristics. Sin ce this model grid i s base d on both mass absorption 
coefficients dKoike et alj|2003l 120061) of several different sam- 
ples and a DHS mod el with a grain siz e of 0.1 fim using the 
optical constants from ISuto et al.l (120061) . we can give neither a 
shape model nor a fixed grain size for the individual models in 
the interpolated grid. 

We performed a minimum x 2 analysis, comparing the inter- 
polated laboratory measurements to the observed forsterite fea- 
tures. The resulting x 2 maps and best-fit curves are plotted in 
Figs|6]and|7] 



From the minimum* 2 analysis we computed confidence in- 
tervals for the iron content and the dust temperature. First, we 
integrated the reduced x 1 distribution with v degrees of freedom 
from x — x 2 /v to infinity. The new distribution, Q(x 2 /v, v), de- 
scribes the probability that a value of x 2 / v or larger is produced 
by random noise. Consequently we searched for a value of * 2 /v 
so that 1 - Q(x 2 /v, v) = 0.997. In Table[3]we list the highest and 
lowest temperature and iron content fractions within the range of 
models below the 3 cr threshold. In the case of AS 205, however, 
the minimum * 2 /v is already larger than the computed 3 <x limit 
as the fit is dominated by systematic errors due to a deviation of 
our band model from the measured profiles at very low temper- 
atures. We therefore cannot give a proper confidence interval for 
the dust temperature or iron content in that system. We discuss 
this result in more detail below. 

Our detailed analysis confirms our initial estimate that the 
observed 69 /mi features are consistent with iron-poor olivine. 
From Fig. |6]it is clear that the olivine dust grains in all sources 
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Table 3. Confidence intervals for iron fraction, grain temperature 
and dust distance to the host star. 



AB Aur 



Star 



Iron fraction [%] Temperature [K] distance [AU] 





mm 


max 


min 


max 


min 


max 


AB Aur 


1.9 


3.5 


74 


273 


16 


221 


HD 100546 


0.1 


0.3 


184 


223 


20 


29 


HD 104237 


0.4 


1.2 


60 


184 


31 


289 


HD 141569 


0.0 


1.2 


107 


>300 


<9 


72 


HD 179218 


0.4 


0.7 


126 


173 


104 


196 


HD 144668 


0.0 


0.4 


130 


224 


25 


74 


IRS 48 


0.1 


0.6 


124 


195 


17 


43 


AS 205 




0.0 




121 


32 





Note: Iron fraction and dust temperature are fitted while the distance of 
the dust to the host stars is estimated based on the temperature (see text) 
The best-fit models are shown in Fig. [7] 



but AB Aur contain at most 1-2% of iron. Though the olivine 
grains in AB Aur can not be iron-free, with a minimum frac- 
tion of 2%, the iron fraction can not be more than 5%. The dust 
temperatures are less well constrained as the iron fraction though 
most fitted temperatures seem to be between 100-200 K. We will 
return to this in the following sections as we discuss the effects 
of temperature and grain size distributions. 

The poorest match to the observed band profile is achieved 
for AS 205. Our best fit model has a band profile which is too 
narrow, and peaks at slightly longer wavelengths compared to 
the observed band, as can be seen from Fig. [7] We believe this 
is a result of the deviation of the linear model of th e FWHM as 
functi on of temperature from the laboratory data by iKoike et all 
(2006J) for low temperatures (<30 K). Lattice distortions in the 
laboratory samples are know n to lead to an in crease in FWHM, 
independent of temperature (llmai et al . 2009). If the grains are 
more distorted than those used in the laboratory experiments a 
higher FWHM would be expected. Also the grain sh ape can in- 
fluenc e the resulting spectrum (see Section 3.1.3; Ima i et all 
2009). If the oli vine grains in AS 2 05 are more spherical than the 
sample used by IKoike et alJ (|2006), they would peak at slightly 
shorter wavelengths (see also Fig. |3). Independent of these un- 
certainties, it is clear that the forsterite grains in AS 205 must 
be extremely iron-poor. Even small fractions of iron would shift 
the forsterite feature to such long wavelengths that effects like 
lattice distortions or grain shape would not explain the observed 
band position any longer. 



4.2. Spectral decomposition: the effect of a temperature 
distribution 

In addition to the analysis of the previous section we also inves- 
tigated the influence of a temperature distribution on the derived 
dust properties. For a more precise analysis of the temperature 
of the forsterite, we fitte d the observed profiles with a weighted 
sum of laboratory data (ISuto etal.ll2006t IKoike et all|2003h at 
different temperatures (50, 100, 150, 200 and 295 KlJ(Fig. [8]). 
The measured flux F^^^m is modeled by 



F 61 - 72m (V = F co M) + K(r dusv A) ■ B A (A, T' dusl ). (2) 



~3 




1 The peak of the emission band in AS 205 is at less than 69 yum and 
therefore can only be fitted with laboratory data in the 8-150 K range 



69 70 
wavelength [|im] 

Fig. 9. AB Aur: The continuum-subtracted forsterite emission at 
69 /im (histogram) over-plotted with the best-fit model (red line) 
which is given by eq. [2] This model is a sum of the following 
components: 50 K 3% iron (red dotted line), 50 K with 4% iron 
(yellow dashed line), 100 K 3% iron (green dash-dotted line) and 
100 K 4% iron (blue dash-three-dots line) and 200 K 4% (purple 
dashed line). Each of the components is multiplied by a fitted 
weighting factor. 



Here F cont is the local continuum (a 2 nd or 3 rd order polynomial), 
and the k values are the mass abs orption coe f ficient s, computed 
from the o ptical constants of the ISuto et all (120061) or directly 
taken from Koike et al. I d2006h pure forsterite sample and B A are 
the Planck function at the dust temperatures T'. .. Finally, the 
w, values are the relative weights of the flux at the different 
temperatures. In addition, we had to specify a dust shape dis- 
tribution in order to calculate the mass absorption coefficients 
from the optical c onstants given bylSuto et a!] (120061) . We choose 
the DHS model dMin et alj|2003l) with a grain size of 0.1 yum 
and a filling factor of / max of 1, which has been found to be 
a good re presentation of ob served silicate profiles in Spitzer 
data (e.g jjuhasz et al.l l2~010). The parameters are fitted with the 
IDL implementat ion of the least squares minimisation rapfit 
(lMarkwardtll2009l) . 

In a first atte mpt we only use the pure forsterite data of 
ISuto et all (120061) in the fit, except for AS 205 where the labora- 
tory data from Koike et all (120031) was taken to account for the 
position of the peak (see Figure[8]for results). All detected bands 
with the exception of AB Aur are very well described under this 
assumption by the model from Eq. [2] The most dominant com- 
ponents for almost all sources are at 1 00-200 K, a result which 
is consistent with the analysis presented in Section |4~T1 The ma- 
jor contribution in AS 205 stems from dust at 20 K. In one case 
(HD 141569) a large spread of temperatures (50-295 K) is found 
which may suggest two separate reservoirs of forsterite in the 
disk. As we show in section l4~4l the shape of the emission band 
in HD 141569 can alternatively be explained by large (-10 jjm) 
grains at a single temperature. The bump in the data at around 
70.4 fim is not a significant feature and consequently we do not 
investigate it further. 

The peak of the emission band in AB Aur is found beyond 
70 fim and can only be fitted if we allow for an admixture of 
some (~ 3-4%) iron. The best fit to the AB Aur spectrum using 
eq. |2]that includes iron is shown in fig. [9] The mass absorption 
coefficients used for the se models were in terpolated from the 
ISuto et all (120061) and fhe lKoike etal] (120031) data. For the inter- 
polation we assumed that the band is shifted linearly with the 
iron content and exhibits the same temperature dependency for 
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Fig. 6. Resulting reduced x 1 distributions from a comparison between our model grid for the 69 //m forsterite feature as a function 
of iron fraction and grain temperature, and the Herschel-observed bands in eight of our targets. Our model grid is an interpolation 
based on several different measurements of both optical constants with the DHS shape model and with absorption coefficients. Thus 
no fixed grain size can be given. Where applicable a grain size of 0. 1 fim was used. 
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Fig. 7. Model fits to the observed 69 /mi forsterite feature. The models correspond to the minimum in the^ 2 distributions shown in 
Fig. [6] Our model grid is an interpolation based on several different measurements of both, optical constants with DHS shape model 
and absorption coefficients. Thus no fixed grain size can be given. Where applicable a grain size of 0.1 fim was used. 
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all different iron fractions up to 10%. Even for AB Aur we find 
that no more than 3-4% iron is needed to obtain a good fit. 

For each source we computed a mass-averaged temperature 
from our fit results. The confidence intervals for the dust tem- 
peratures were derived from a x 2 analysis similar to the proce- 
dure we discussed in section |4~T1 Since we only have laboratory 
data at five different temperatures, and the effects of nonzero iron 
fractions < 1 % cannot properly be taken into account as all labo- 
ratory measurements were made using pure forsterite, the values 
derived here may suffer from larger systematic effects than the 
ones given in Section 14.11 For AS 205 and AB Aur we cannot 
compute confidence intervals as even the best-fit value already 
has too large of a reduced x 2 - in the case of AS 205, the mis- 
match between the observed band profile and the best-fit model 
can be explained as a consequence of insufficient laboratory data 
concerning the optical properties at temperatures of <50 K. For 
AB Aur we have to manually select the iron content in the data 
used to fit the temperature distribution. Since the iron content is 
not a parameter in the fit in this model (see eq. |2]l the minimal 
X 1 achieved in the fit is too large for the computation of a 3 cr 
confidence interval. 

We summarise our fit results in Table |4] All detected 69 /mi 
bands except the ones in AS 205 and AB Aur can be fitted with 
pure or iron-poor and warm (150-200 K) forsterite. In AS 205 
the most dominant component is found at -20 K, while AB Aur 
requires some percent of iron at a temperature of ~50 K. 

The temperatures, and therefore the radial distance estimates 
based on them, are consistent with those derived from the x 1 
maps in section 14.11 Only AS 205 shows strong deviations be- 
tween the results of the two different methods. This is most likely 
attributed to the poor quality of the fit in either of the two sce- 
narios. In fact, the minimum^- 2 in the fit to the observed profile 
in AS 205 is too large to derive a confidence interval with the 
method described in section |4~T1 

4.3. Distance of the forsterite grains from their host stars 

With the temperature of the forsterite grains derived from the 
profiles of the 69 //m emission band (see Tab. [3] and Tab. |4| we 
can estimate the distance of the dust particles from their host 
star. For this estimate we have to assume that the dust grains 
are located in an environment that is optically thin, in both opti- 
cal wavelengths where the star emits and IR wavelengths where 
the dust grains re-emit, an assumption that is very likely invalid 
for almost all sources in our sample. In that case the estimated 
distances must be interpreted as upper limits. More reliable dis- 
tances would require full radiative transfer modeling for each of 
the sources, which is beyond the scope of this paper. 

The equilibrium temperature of a dust grain in an optically 
thin environment depends on the optical properties not only 
at the wavelengths where it emits but also at the wavelengths 
around the peak of the stellar emission, where it receives most of 
the radiation. Since the forsterite grains in protoplanetary disks 
are likely part of larger, porous aggregates mixed together with 
other materials, these properties will be very different from those 
of pure forsterite grains. However, the exact composition and 
structure of these dust aggregates cannot be determined from the 
infrared spectra in the Spitzer IRS and Herschel PACS wave- 
length ranges since most of the possible constituents dominating 
the opacities at optical wavelengths ( e.g. carbon, iron) have no 
emission bands in this spectral region. iLisseet all (120091) present 
an approach, where it is assumed that the dust aggregates in 
circumsteilar disks have optical properties similar to those in 
cometary material in the solar system for which we know the 



temperature-dis tance relation. We therefore use the relation from 
iLisse et aD (I2009L section 2.2.6) and solve for the distance of the 
forsterite particles from their respective host star. 

We estimate the location of the forsterite grains in the disks 
based on the dust temperatures derived i n sections!4. 1 I and l4.2l us- 
ing the temperature-distance relation bv lLisse et al.l (12009). The 
derived distances are listed in Tables [3] and [4] The distance esti- 
mates based on the two different methods used to derive the dust 
temperature agree reasonably for most of the sources with min- 
imum distances of 10-30 AU. The larger values in HD 179218 
are easily explained by the larger luminosity of the star com- 
pared to the rest of our sample. In HD 104237 and AB Aur, 
compositional effects play a more significant role, thus explain- 
ing the different results for these sources in Tables |3]and |4] The 
emission profile of the 69 /mi band in AS 205 is not well-fitted 
in either of our two approaches, leading to a large difference in 
the computed distance of the forsterite from the central star, de- 
pending on the method used to determine it. 

The most notable difference between the results of the two 
methods we used in Sections [4.11 and l4~2l to derive the tempera- 
ture is the size of the confidence interval. The spectral decompo- 
sition method in Section 14721 generally produces much narrower 
confidence intervals on the temperature and thus the radial dis- 
tribution of the forsterite. The reason for this is that in the spec- 
tral decomposition we only use pure forsterite and even a small 
amount of iron could significantly decrease the temperature of 
the best-fit solution (see section lXTT i. The effect is also illustrated 
in Figure [6] where the x 2 contours extend to lower temperatures 
as we move to higher iron fractions. 

In fact all the distance values given in Tables[3]and|4]are es- 
timates and should be taken with caution. For the computation 
we have assumed that the forsterite is located in an optically thin 
environment. This assumption may be violated in protoplane- 
tary disks li ke the ones observed in our sample. As shown for 
example by M ulders et al.l (1201 ll) the temperature profile inside 
an optically thick disk is dramatically different from that of the 
optically thin disk atmosphere. 

A more precise determination of the location of the forsterite 
in the protoplanetary disks requires detailed radiative transfer 
modeling which is beyond the scope of this paper. Such mod- 
eling would also greatly benefit from a more complete set of 
optical constants, covering olivines with low iron content at dif- 
ferent temperatures. Radiative transfer modeling is very likely 
to decrease the derived distance of the forsterite reservoirs from 
their central stars, unless the grain s are located at the i nner edge 
of a disk (as in HD 100546, see iMulders et al1l201 lh . The dis- 
tances derived in this paper can be regarded as upper limits. 

4.4. Grain size effects 

The effect of grain size of the emitting dust on the position and 
shape of the 69 /mi band has been discussed in Section |3~T1 and 
illustrated in Figure [4] The largest grains that produce a 69 /mi 
band with a profile that is compatible to our detections have radii 
of about 10 /an. 

If we take large grains (1-10 /mi) into account, the tempera- 
ture estimate we derive from the peak position and FWHM (sim- 
ilar to Figure|5]l is slightly lower than for small grains. Assuming 
that almost all the dust consists of grains with radii of 10 /mi 
allows a reasonable fit to all detected bands with only one tem- 
perature component at around 100 K. However, for most of the 
objects the fit is significantly worse than the one presented in 
Figure [8] where we used smaller grains. 
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Table 4. Fit results of the spectral decomposition of the detected forsterite bands. 









Relative mass fractions 




M.A.T 


distance 


Star 


@ 50 K 


@ 100 K 


@ 150 K 


@ 200 K 


@ 295 K 


[K] 


AU 


HD 100546 


0.00 


0.11 


0.22 


0.63 


0.04 


166-204 


24- 36 


IRS 48 


0.00 


0.00 


0.45 


0.55 


0.00 


156-189 


19- 27 


HD 144668 


0.00 


0.00 


0.93 


0.07 


0.00 


135-183 


37- 69 


HD 179218 


0.00 


0.00 


U.JZ 


U.Oo 


0.00 


174 - 186 


90 - 103 


HD 141569 


0.00 


0.00 


0.37 


0.23 


0.40 


148 - 248 


13- 38 


HD 104237 


0.00 


0.00 


0.00 


1.00 


0.00 


180-211 


23- 32 




@ 8.K 


@ 20.K 


@ 50.K 


@ 100.K 


@ 150.K 






AS 205 


0.00 


0.78 


0.22 


0.00 


0.00 


~ 30 


~ 516 




@ 50 K, 3% Fe 


@ 50 K, 4% Fe 


@ 100 K, 3% Fe 


@ 100 K, 4% Fe 


@ 200 K, 3% Fe 






AB Aur 


0.75 


0.00 


0.10 


0.15 


0.00 


~ 63 


~ 304 



Note: Listed are the relative mass fractions of the forsterite components at different temperature or with different iron mass fractions. Also listed 
is the confidence intervals for the mass-average-temperature (M.A.T.) of the fitted dust components. The distance intervals are computed based on 
the ranges of dust temperature (see text). 



HD 141569 




68 69 70 
wavelength [\im] 

Fig. 10. The 69 pm band in HD 141569 fitted with a combination 
of the emission from different grain sizes at 150 K. Only the 
emission from grains of 10 pm diameter contributes to the best- 
fit model. 



A more detailed analysis by fitting a weighted sum of dif- 
ferent grain sizes to the data (sim ilar to Eq. |2} is done for each 
of the temperatures at which lSuto et al.l (120061) published optical 
constants (50, 100, 150, 200 and 295 K). The best results are typ- 
ically found for 100-150 K dust with no significant contributions 
(mass fraction < 10~ 3 ) from 5-8 pm grains. AB Aur cannot be 
fitted with pure forsterite and we have no laboratory data avail- 
able from which we could derive the cumulative effect of iron 
and grain size. The band in AS 205 is at very short wavelength 
and very narrow. None of the fits with large grains could describe 
the observed emission band. Only for HD 141569 do we find in- 
dications for a significant contribution from 10 pm grains. For 
this system the 69 pm band can be modeled with emission from 
10 pm grains only (Fig. ITOb at a temperature of 150 K without 
any significant change in fit quality compared to the model we 
discussed in the previous subsection. 

4.5. Comparison to Spitzer IRS data 

Most of our targets were observed with the InfraRed 
Spectrograph (IRS) onboard the Spitzer Space Telescope. To 



compare the results of our analysis of the Herschel spectra with 
those of the forsterite emission bands s een in the Spitzer o bser- 
vations, we use the s pectra published bv lJuhasz et alj d2010l) and 
Olofs son et al.l d2009h . For the targets where we detected a 69 pm 
band, and those with marginal Spitzer detections, we carefully 
re-reduce d the Spitzer spectra , using the same data reduction ap- 
proach as Juha sz et"ai1 (120 1 Ot) . but with the latest calibration ver- 
sion (S 18. 18.0). Our re- processed spectra on ly m arginally differ 
from t hose reported by Juh asz et"aTI (|2010) and Olofsso n et al.l 



(120091) . However, we report the detection of forsterite bands be- 
tween ~19 t o 33 pm in the Spitz er spectrum of HD 141569, not 
reported by Juhasz ^tall (l2010l). Interestingly, this source was 
classified by iMeeus et al.l (1200 ll) as a Group lb source, i.e. a 
strongly flaring disk with no silicate emission bands. 

We searched for correlations between the strength of the 
forsterite bands seen in the Spitzer data and the 69 pm band. For 
this, we analysed the forsterite features seen in the Spitzer spec- 
tra in an identical fashion to our analysis of the 69 pm band. For 
all the forsterite bands we computed the ratio of the Lorentzian 
and the continuum polynomial at the position of the peak of the 
best-fit Lorentzian. The comparisons between the 16 pm and 
33 pm peak over continuum values derived from the Spitzer data 
and the 69 pm peak over continuum values are shown in Figures 
Q~T|and[T2] We compare the Herschel results with the 16 and 
33 pm bands as those have very little confusion with other emis- 
sion bands, and probe the warmest and coldest dust detectable in 
the Spitzer wavelength range, respectively. 

The comparison between the 16 pm and 69 pm bands does 
not show any obvious correlation, implying that there is no 
straightforward observational correlation between the emission 
from the warmest and the coldest forsterite grains as probed by 
the 16 pm and, respectively, the 69 pm band. There appears to 
be a slightly stronger correlation with the 33 pm band: those 
sources which show the strongest 69 pm band also show an on 
average stronger 33 pm band. A correlation between the 33 pm 
and the 69 pm feature might be expected due to the relative warm 
temperatures of the forst erite grains (~ 150- 200 K, see Table 0). 
However, as discussed bv lSturm et~ail(l2010l) for HD 100546, the 
33 pm over 69 pm band ratio is consistent with a ~70 K grain 
temperature, rather than the ~ 200 K deriv ed from the 69 fea - 
ture. The solution to this was proposed bv lMulders et al.l (1201 ll) . 
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Fig. 11. Comparison of the peak/continuum of the 69yum to the 
\6pm bands in our sample. The sources with a firm detection of 
the 69 pm band are labeled. 
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Fig. 12. Comparison of the peak/continuum of the 69 pm to the 
33 /mi bands. Sources with a 69 pm detection are labeled. 



who showed that the 33 pm band can be substantially suppressed 
relative to the 69 pm feature by optical depth effects. These 
effects could blur any straightforward correlation between the 
Spitzer and Herschel observations. 

To check whether any of the other disks could have similar 
optical depth effects as those seen in in the disk of HD 100546, 
we plotted in Fig.[T3]the observed 33 pm over 69 pm band ratio 
against the expected ratios based on the temperatures listed in 
Table [4] Except for AS 205, all disks must be optically thick 
at mid-IR wavelengths to suppress the forsterite bands in the 
Spitzer IRS wavelength range. For these disks, the emission seen 
with Spitzer arises from the disk surface, while the bulk of the 
69 pm emission originates from the disk interior. Most of the 
emission should come from a disk region where dust tempera- 
tures in the range of 150-200 K are reached, which for a typi- 
cal Herbig Ae system is ~10 AU from the central star. Only for 
AS 205 are our findings are consistent with a complete lack of 
warm forsterite grains. In that case, the emission originates from 



Fig. 13. Comparison between the measured 33 pm over 69pm 
peak fluxes in the forsterite bands and the expected ratio's based 
on the derived mass averaged temperatures of the forsterite 
grains as listed in Table [4] The dashed line is the expected be- 
haviour if no optical depth effects play a role. 
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Fig. 14. Spectral type versus signal-to-noise ratio. Sources with 
a detected 69 pm forsterite band are shown as red squares. The 
S/N is computed as the average flux divided by the mean point- 
to-point difference in the spectra. 



a reservoir of cold grains with very low iron content located ei- 
ther in a region far out in the disk, or in the cold interior, which 
is shielded from direct stellar radiation. Note that as this system 
is a T Tauri star, low temperatures are reached closer to the cen- 
tral star than in the Herbig Ae systems. The physical location of 
the forsterite grains might therefore not be that different from the 
~10 AU we expect for the intermediate-mass systems. 

In contrast to the amorphous silicate grains, no evidence 
for substantial growth o f the forsterite cry stals is found from 
Spitzer observations (e.g. lJuhasz et al.l l2010). The observed mid- 
IR emission bands are all consistent with emission from small 
(<1 pm) forsterite grains. Also in our Herschel observations no 
conclusive evidence is found for the presence of larger grains. 
The only source where indications of larger (10 pm) grains are 
found is HD 141569, though the 69 pm band can also be mod- 
eled with small grains only. As the far-IR observations tend to 
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probe closer to the disk midplane, if grain growth and grav- 
itational settling of the larger grains occur in these disks, our 
Herschel observations should be more sensitive to a population 
of larger grains. As we do not detect these larger grains, these re- 
sults are consistent with the notion that the small forsterite grains 
are embedded into larger aggregates consisting o f mainly amor- 
phous silicates or carbonaceous material. As iMin et al.l ([2008) 
show for aggregates, materials that have a very low abundance 
appear spectroscopically as if they were in very small grains, 
while more abundant materials appear spectroscopically to re- 
side in larger grains. 

We also investigated whether a correlation exists between 
the presence of th e 69 fim band and the global disk shape. 
iMeeus et alJ d200lh divided the Herbig Ae/Be systems into 2 
groups: Group I, having a large far-IR excess corresponding with 
a strongly flaring disk, and group II, having a SED consistent 
with a flattened disk structure. HD 104237 is the only group II 
source with a forsterite detection; all other sources are classi- 
fied as group I. This might, however, be an observational bias. 
Due to the SED shape of the group II sources, the far-IR fluxes 
are systematically lower than those of the group I sources and, 
therefore, the S/N of the spectra are typically lower. An exam- 
ple of this is HD 150193, a source with one of the strongest 
forsterite features at mid-IR wavelengths, but with no detected 
69 fim band. Due to the low far-IR flux of this source, the up- 
per limit we can place on the forsterite band would still be con- 
sistent with a relatively strong feature, like the one observed in 
HD 179218. 

4.6. Correlation of detection rate with S/N and spectral type 

We carefully looked for a correlation of the detection of a 69 fim 
band with the S/N. To account for the complex continua in our 
sample, we computed the S/N as the average flux divided by 
the mean point-to-point difference in the spectra. Among the 16 
sources with the highest S/N, we find five detections, compared 
to three detections among the 16 sources with lower S/N. This 
indicates a weak correlation of the detection rate with S/N. 

Remarkably, almost all stars with a detection of the 69 fim 
forsterite feature have spectral types of B9-A9 (see fig. [14} . The 
only exception is AS 205, which has a spectral type of K5. In to- 
tal we have detected the 69 fim band in 7 out of 21 objects with 
spectral types earlier than F compared to only 1 detection among 
1 1 objects of spectral types F and later. This seems to indicate 
a correlation of the detection rate with spectral type. However, 
the correlation is not significant. We show the spectral type dis- 
tribution in our sample and the distribution of the detections of 
the 69 fim forsterite band in Fig. Q3] The number of bins in the 
histog ram was chosen following the strategy presented in Hogg 
(2008). Assuming the 8 detections of the 69 fim forsterite band 
were randomly distributed among the 32 sources in our sample, 
the probability to find the observed distribution would be -8%. 

We also carefully checked if there are significant deviations 
in the S/N distribution within our target sample for stars with 
spectral types earlier or later than F0. A non-parametric Mann- 
Whitney two-sample U-test indicates that the probability that 
the S/N distribution is the same for both early- and late-type 
stars is P — 0.38, or formulated differently, the probability that 
our Herbig Ae/Be sample has a different S/N distribution than 
the T Tauri sample has only 1 cr significance. Thus, we have 
not detected a significant S/N bias favouring either early or late 
spectral types in our sample. 

AS 205 may show up as the only 69 /jm detection in a spec- 
tral type later than A because of the geometry and orientation 



25 




spectral type 

Fig. 15. The spectral type distribution in the observed sample 
(blue) and the distribution of the 69 fim detections (red). 

of this particular disk. lAndrews et al.l (120091) find that the scale 
height of the protoplanetary disk in AS 205 is the largest in 
their sample, while the mass is at the lower end of the spectrum. 
This implies that the dust density is low compared to the other 
sources, allowing for an optically thin medium and therefore en- 
abling us to detect all of the forsterite mass at 69 fim while in 
more massive/dense disks only a fraction would be visible. This 
is supported by our finding that we do not need to take opti- 
cal depth effects into account when comparing the Herschel and 
Spitzer data (see Fig.[T3l). Also, the large scale height indicates a 
strongly flaring disk, which will intercept a higher portion of the 
central star's radiation compared to flatter disks. Thus, the disk 
of AS 205 might be warmer than those of stars with comparable 
spectral types. Finally, the disk of AS 205 is seen nearly face-on 
( lAndrews et al.ll2009l) . maximising the projected surface which 
is seen by Herschel. 

5. Conclusions 

We searched for the 69 fim forsterite emission band in 32 disk 
sources (see Table[TJ and identified it in only eight of the objects. 
Several more spectra show some emission features between 69 
and 70 fim which can not clearly be attributed to forsterite. This 
detection rate (even including the features with uncertain iden- 
tification) is much lower than the one for the forsterite bands in 
the Spitzer IRS wavelength range of the same sample. 

This work provides the first overview of the presence and 
properties of the 69 fim band in a sample of disks that covers 
various ages, disk masses and stellar properties. Only one detec- 
tion of the 69 pm band in a di sk source was repor ted from an 
ISO observation (HD 100546, iMalfait et alj|1998l) . Thanks to 
the unprecedented sensitivity and spectral resolution of Herschel 
PACS we can for the first time derive the temperature and com- 
position of the crystals from the shape and position of the emis- 
sion band. 

5.1. Iron content, size, temperature and location of the 
forsterite grains 

We compared the observed 69 fim forsterite emission band with 
laboratory measurements of this feature for different iron frac- 
tion and grain temperature. In accordance with the ISO detection 
of the 69 fim band in HD 100546 we find that all but one of the 
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detected bands in our sample are consistent with emission from 
iron-poor grains. For one source, AB Aur, a 3—4% iron fraction 
is needed in order to explain the position of the detected emis- 
sion band. In all other objects pure forsterite can account for the 
69 pm features although a small amount of iron (< 1 %) cannot be 
ruled out. Our results are consistent with the recent detection of 
the 69 pm band in the debris disk system j3 Pictoris, which shape 
and positi on is consistent with olivine with an iron fraction of at 
most 1 % dde Vries et al.ll2012l) . 

Our analysis shows no conclusive evidence for larger (5- 
10 pm) forsterite grains. The observed 69 pm bands can be mod- 
eled with emission from small (1 /mi or less) forsterite grains. 
These inferred small grain sizes are consistent with the grain 
sizes for the forsterite cr ystals derived from Spitzer observations 
(e.g. lJuhasz et al.l l2010). The only source where indications of 
larger (10 pm) grains are found is HD 141569, though the ob- 
served 69 pm feature can also be modeled with small grains only. 

From the analysis of the position and shape of the emission 
bands we derive a best-fit grain temperature for pure forsterite 
grains of around 150-200 K for all sources except for AB Aur 
and AS 205. The latter source requires very cold dust grains at 
temperatures of 30 K or less to explain the observed band posi- 
tion, though an exact temperature could not be determined due 
to the sparsity of laboratory measurements of the optical prop- 
erties of forsterite at these very low temperatures. In the case 
of AB Aur, no pure forsterite can explain the observed spectral 
feature. Determining grain temperatures using iron-containing 
olivine grains leads to similar results though the derived temper- 
atures are less well-constrained due to systematic uncertainties 
(see Section 3.1) caused by the interchangeability of the effects 
of grain temperature and iron content on the band position and 
shape. 

We used the derived grain temperatures to estimate the radial 
distance from the forsterite grains to their respective host stars. 
We find that the forsterite grains we observe with Herschel are 
located at least 10 AU or farther from the central star. In case 
of pure forsterite grains, we find the forsterite for most of our 
detections (except AB Aur and AS 205) to be located in a rel- 
atively narrow ring in the circumstellar disk. Note that for our 
simple distance estimation we assume an optically thin environ- 
ment and that all of the values will be modified when taking 
radiative transfer effects in the disks with wavelength-dependent 
optical depth into account. The distances estimates in this pa- 
per should, therefore, be regarded as upper limits. In the case of 
HD 100546, an analysis of the forsterite grain lo cation using a 
full ra diative transfer model has been published (tMulders et alJ 
l201lh . The derived location in a narrow ring of around 13-20 AU 
radius by|Mulders "eTaII(l20T ]]) is at least qualitatively in agree- 
ment with the narrow location we derive of ~20-30 AU. 

5.2. Formation history of the forsterite grains 

With the constraints on the grain temperature and composition 
we derived from the Herschel PACS spectra, we can specu- 
late on formation history of the forsterite grains. The very low 
iron content in the crystalline silicates supports a crystal forma- 
tion scenario through co ndensation from th e gas phase at high 
(-1500 K) temperatures dGail & Hoppdl2010l) . This is in appar- 
ent contradiction to the relatively low grain temperatures we de- 
rive for the forsterite, as was already noted from the analysis 
of ISO and Spitzer spectra (e.g lJuhasz et al.ll2010t iMalfait et alJ 
1998). Several scenarios have been proposed to explain this dis- 
crepancy, such as radial mixing from the hot inner d isk to the 
outer disk regions (e.g. Bockelee -Morvan et al.l |2000T) or tran- 



sient heating events like shocks ( Harker & Deschl [20021) tem- 
porarily heating up circumstellar material in disk regio ns other- 
wise too cold for forsterite formation to occur (see also lHenningj 
120101 for a review on dust formation and processing). 

Given the large uncertainties in the temperature estimate and 
the fact that our assumption of an optically thin medium may 
be violated we can only speculate whether any radial mixing 
could have occurred. Since the temperature gradient in optically 
thick disk s is much steeper t han in an optically thin environ- 
ment (e.g. Mul ders et al.ll2.01 lb . not only will the best estimate 
for the dust distance become smaller but also the confidence in- 
terval narrower if we include radiative transfer effects. Thus it 
seems likely that the forsterite we detected is mostly located in 
narrow rings/tori. The distribution produced by large-scale ra- 
dial mixing, however, i s expected to be muc h more radially ex- 
tended (e.g. lGaiill200ll) . lJuhasz et al.l (120101) also argue against 
radial mixing due to a lack of crystalline dust species other than 
forsterite at lower temperatures, which should have higher mass 
fractions than forsteri te at the ou ter disk regions, according to 
radial mixing models (Gail 2Q0lb. 

For one source in our sample, HD 100546, we seem to be 
able to rule out the radial-mi xing scenario. For t his system a 
full radiative transfer model (M ulders et al.ll20TTb is available 
which shows that the forsterite is concentrated in a rather nar- 
row ring. Since the temperatures in most of the sources are sim- 
ilar to HD 100546, it seems reasonable to assume that the dis- 
tribution of the forsterite grains in those sources is also similar 
to that of HD 100546. Thus it seems unlikely that radial mix- 
ing can explain the narrow distribution of forsterite found in 
in the disks of our sample. As an alternative to radial mixing 
iBouwman et al.l (120031) propose a connection between the forma- 
tion of the crystalline silicates on site and the creation of a gap 
in the disk of HD 100546 through the gravitational interaction 
with a massive planet at around 10 AU. A planet opening a gap 
in the disk, would prod uce shock waves through tidal interac- 
tion w ith the disk (e.g. lLin & Papaloizou 1980; Boss & Durisenl 
2005), which might be enough to heat the circumstellar ma- 
terial Jo_Jiighenough temperatures for crystallisation to oc- 
cur dDesch et al.ll2005h . Planetary companions are not the only 
possible cause for Shockwaves in protoplanet ary disks, which 
could also be induced by disk instabilities (e.g. Harker & Desch 
2002) . Given that all of our targets are still gas-ric h (Meeus et al 
2012; lFukagawaet~al1l2010l; iBrown et alJl2012t lAndrews et al. 
2009, and references therein), small dust grains should be effi- 
ciently heated to high temperatures at disk radii of up to at least 
-10 AU. 

Alternatively, collisions of large, possibly differentiated ob- 
jects (plan e tesima ls), could also produce crystalline olivines. 
iLisse et ail (120091) found that the Spitzer spectrum of the de- 
bris disk around HD 172555 is dominated by silica emission 
and shows possible features arising from SiO gas, which the 
authors explained as evidence for hypervelocity collisions of 
planetesimals. However, contrary to the system of HD 172555, 
only v ery small amounts of silica are reported by Juhas z~et~al1 
(12010 ) for our targets where we detect the 69 pm forsterite emis- 
sion band. Also no clear indication of SiO gas can be found 
in the residuals of the IRS spectra. Further, such c ollisions are 



less li kely in the gas disks found in our samp le (Meeus et al 
2012; Fukaga wa et al] 120 lot iBrown et all 120 121; Andrews et al. 
2009, and references therein) than in a gas-poor debris disk 
like HD 172555 since the interaction between disk and plan- 
etesimals will dram atically reduce the eccentricity of orbits 
dBitsch & Klevll2010l) . 
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Furthermore, the formation of olivines in larger, differenti- 
ated objects is a scenari o that would result in a much higher iron 
content than we found dGail & Hoppdl2010l) . It is interesting to 
note that only for a few debris disk sources is there solid evi- 
dence for olivine with a high (~20%) iron fraction, suggesting 
that the olivine material in these debris disks o riginated from 
larger differentiated bodies (Olofs son et 51 12012). This is also 
consistent with the recent detection of the 69 //m band in the 
young debris disk system /3 Pictoris, which shape and position 
is consistent with ol ivine with an iron fraction of at most 1 % 
dde Vries et al.ll20"l2l) . It seems that the bulk of the olivine dust 
in protoplanetary and young debris disks is formed through an 
equilibrium condensation process at high temperatures, and that 
only during the later debris disk phases a substantial amount of 
iron-rich crystalline silicates is produced through disruptive col- 
lisions of differentiated bodies. 

From the evidence at hand we conclude that it is unlikely, al- 
though not impossible, that the crystalline silicates in our sample 
to be produced by large hypervelocity collisions of differentiated 
objects. 

5.3. Correlation with mid-IR and disk- or host star properties 

No correlation between the detection or strength of the 69 /mi 
feature and the 16 /mi bands was found. A possible weak cor- 
relation between the 69 /mi and the 33 /mi bands could exist: 
sources with the strongest 69 //m band tend on average to be 
those sources with the strongest 33 //m bands. Note, however, 
that a detection at 69 //m does not automatically imply the pres- 
ence of forsterite bands at the shorter wavelengths. Examples 
for this are AB Aur, IRS 48 and AS 205 which have no de- 
tectable forsterite bands at mid-IR wavelengths and still show- 
ing a 69 yum band. 

No conclusive correlation between the existence of the 
69 /mi band and the disk geometry has been found. Although 
only one of the eight forsterite detections was found in a flattened 
disk (HD 104237), the possibility exists that this is an observa- 
tional bias due to the systematically lower fluxes of sources with 
flattened disks compared to those with strongly flaring disks. 

In Figures [14] and [15] we show that the detection rate of the 
forsterite 69 //m band is highest among early spectral types with 
high S/N. The correlation between spectral type and detection 
rate indicated by this finding is not significant. The chance of 
finding the observed distribution of detections in our sample by 
chance is ~8%. The disk of AS 205 has geometric and compo- 
sitional properties such as large scale height, low density, and 
near face-on orientation (lAndrews et al.ll2009l) . unusual among 
the disks of T Tauri stars. These properties increase the chance 
that Herschel would be able to detect the 69 /mi band. Therefore, 
we conclude that the indicated correlation of detection rate and 
spectral type is most likely an observational bias. 
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Fig. 8. Results of fitting eq.[2]with pure forsterite to the detected 69 /im forsterite emission bands. The measured spectrum is plotted 
as a black line, over-plotted with the best-fit model (continuous red line). The contribution of crystals at different temperatures is 
shown as follows: The dotted red line is the 50 K component, dashed yellow represents 100 K, dash-dotted green stands for 150 
while the dashed-thr ee-dot blue curve represents the 200 K sample and dashed purple the 295 K data. In all fits except for AS 205 
we use the data from lSuto et al.1 ([2006) with the DHS shape model an d a grain size of .1 fjm. In AS 205 the position of the peak 
indicates very low temperatures so we used the laboratory data from dKoike et alJl2006h which include 8- (red, dotted) and 20 K 
(dashed, yellow) samples. The profile of AB Aur cannot be fitted with pure forsterite. In this case the admixture of some percent 
iron greatly improves the fit (see fig. |9) . 
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Fig. A.l. Comparison of the unmodified (black) and the noise-filtered (red) spectrum in the left column. The middle column shows 
the noise-filtered spectrum (black) overplotted with the best-fit model (red). The model is a sum of a 3 rd order polynomial and a 
Lorentzian. In the right column we show a comparison of the residuals to the Lorentzian from the model. The error bars on the 
residuals are corresponding to the standard deviation of the residuals. In all sources shown in this figure the 69 fim forsterite band 
has been detected. 
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Fig. A.2. Same as fig. IA. 1 1 In all sources in this figure the 69 //m forsterite band has been detected. 
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Fig. A.3. Same as fig. IA. 1 1 In the sources in this figure, the 69 jum forsterite band has not been detected. 
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Fig. A.4. Same as fig. IA. 1 1 In the sources in this figure, the 69 jum forsterite band has not been detected. 
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Fig. A.5. Same as fig. IA. 1 1 In the sources in this figure, the 69 jum forsterite band has not been detected. 
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Fig. A.6. Same as fig. IA. 1 1 In the sources in this figure, the 69 jum forsterite band has not been detected. 
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Fig. A.7. Same as fig. IA. 1 1 In the sources in this figure, the 69 jum forsterite band has not been detected. 
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Fig. A.8. Same as fig. IA. 1 1 In the sources in this figure, the 69 jum forsterite band has not been detected. 



